3.1662 \(\int \frac{1}{(d+e x)^{5/2} (a^2+2 a b x+b^2 x^2)^2} \, dx\)

Optimal. Leaf size=200 \[ \frac{105 b^{3/2} e^3 \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 (b d-a e)^{11/2}}-\frac{105 b e^3}{8 \sqrt{d+e x} (b d-a e)^5}-\frac{35 e^3}{8 (d+e x)^{3/2} (b d-a e)^4}-\frac{21 e^2}{8 (a+b x) (d+e x)^{3/2} (b d-a e)^3}+\frac{3 e}{4 (a+b x)^2 (d+e x)^{3/2} (b d-a e)^2}-\frac{1}{3 (a+b x)^3 (d+e x)^{3/2} (b d-a e)} \]

[Out]

(-35*e^3)/(8*(b*d - a*e)^4*(d + e*x)^(3/2)) - 1/(3*(b*d - a*e)*(a + b*x)^3*(d + e*x)^(3/2)) + (3*e)/(4*(b*d -
a*e)^2*(a + b*x)^2*(d + e*x)^(3/2)) - (21*e^2)/(8*(b*d - a*e)^3*(a + b*x)*(d + e*x)^(3/2)) - (105*b*e^3)/(8*(b
*d - a*e)^5*Sqrt[d + e*x]) + (105*b^(3/2)*e^3*ArcTanh[(Sqrt[b]*Sqrt[d + e*x])/Sqrt[b*d - a*e]])/(8*(b*d - a*e)
^(11/2))

________________________________________________________________________________________

Rubi [A]  time = 0.140279, antiderivative size = 200, normalized size of antiderivative = 1., number of steps used = 8, number of rules used = 4, integrand size = 28, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.143, Rules used = {27, 51, 63, 208} \[ \frac{105 b^{3/2} e^3 \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 (b d-a e)^{11/2}}-\frac{105 b e^3}{8 \sqrt{d+e x} (b d-a e)^5}-\frac{35 e^3}{8 (d+e x)^{3/2} (b d-a e)^4}-\frac{21 e^2}{8 (a+b x) (d+e x)^{3/2} (b d-a e)^3}+\frac{3 e}{4 (a+b x)^2 (d+e x)^{3/2} (b d-a e)^2}-\frac{1}{3 (a+b x)^3 (d+e x)^{3/2} (b d-a e)} \]

Antiderivative was successfully verified.

[In]

Int[1/((d + e*x)^(5/2)*(a^2 + 2*a*b*x + b^2*x^2)^2),x]

[Out]

(-35*e^3)/(8*(b*d - a*e)^4*(d + e*x)^(3/2)) - 1/(3*(b*d - a*e)*(a + b*x)^3*(d + e*x)^(3/2)) + (3*e)/(4*(b*d -
a*e)^2*(a + b*x)^2*(d + e*x)^(3/2)) - (21*e^2)/(8*(b*d - a*e)^3*(a + b*x)*(d + e*x)^(3/2)) - (105*b*e^3)/(8*(b
*d - a*e)^5*Sqrt[d + e*x]) + (105*b^(3/2)*e^3*ArcTanh[(Sqrt[b]*Sqrt[d + e*x])/Sqrt[b*d - a*e]])/(8*(b*d - a*e)
^(11/2))

Rule 27

Int[(u_.)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[u*Cancel[(b/2 + c*x)^(2*p)/c^p], x] /; Fr
eeQ[{a, b, c}, x] && EqQ[b^2 - 4*a*c, 0] && IntegerQ[p]

Rule 51

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)*(c + d*x)^(n + 1
))/((b*c - a*d)*(m + 1)), x] - Dist[(d*(m + n + 2))/((b*c - a*d)*(m + 1)), Int[(a + b*x)^(m + 1)*(c + d*x)^n,
x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && LtQ[m, -1] &&  !(LtQ[n, -1] && (EqQ[a, 0] || (NeQ[
c, 0] && LtQ[m - n, 0] && IntegerQ[n]))) && IntLinearQ[a, b, c, d, m, n, x]

Rule 63

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - (a*d)/b + (d*x^p)/b)^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rubi steps

\begin{align*} \int \frac{1}{(d+e x)^{5/2} \left (a^2+2 a b x+b^2 x^2\right )^2} \, dx &=\int \frac{1}{(a+b x)^4 (d+e x)^{5/2}} \, dx\\ &=-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}-\frac{(3 e) \int \frac{1}{(a+b x)^3 (d+e x)^{5/2}} \, dx}{2 (b d-a e)}\\ &=-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}+\frac{\left (21 e^2\right ) \int \frac{1}{(a+b x)^2 (d+e x)^{5/2}} \, dx}{8 (b d-a e)^2}\\ &=-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}-\frac{21 e^2}{8 (b d-a e)^3 (a+b x) (d+e x)^{3/2}}-\frac{\left (105 e^3\right ) \int \frac{1}{(a+b x) (d+e x)^{5/2}} \, dx}{16 (b d-a e)^3}\\ &=-\frac{35 e^3}{8 (b d-a e)^4 (d+e x)^{3/2}}-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}-\frac{21 e^2}{8 (b d-a e)^3 (a+b x) (d+e x)^{3/2}}-\frac{\left (105 b e^3\right ) \int \frac{1}{(a+b x) (d+e x)^{3/2}} \, dx}{16 (b d-a e)^4}\\ &=-\frac{35 e^3}{8 (b d-a e)^4 (d+e x)^{3/2}}-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}-\frac{21 e^2}{8 (b d-a e)^3 (a+b x) (d+e x)^{3/2}}-\frac{105 b e^3}{8 (b d-a e)^5 \sqrt{d+e x}}-\frac{\left (105 b^2 e^3\right ) \int \frac{1}{(a+b x) \sqrt{d+e x}} \, dx}{16 (b d-a e)^5}\\ &=-\frac{35 e^3}{8 (b d-a e)^4 (d+e x)^{3/2}}-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}-\frac{21 e^2}{8 (b d-a e)^3 (a+b x) (d+e x)^{3/2}}-\frac{105 b e^3}{8 (b d-a e)^5 \sqrt{d+e x}}-\frac{\left (105 b^2 e^2\right ) \operatorname{Subst}\left (\int \frac{1}{a-\frac{b d}{e}+\frac{b x^2}{e}} \, dx,x,\sqrt{d+e x}\right )}{8 (b d-a e)^5}\\ &=-\frac{35 e^3}{8 (b d-a e)^4 (d+e x)^{3/2}}-\frac{1}{3 (b d-a e) (a+b x)^3 (d+e x)^{3/2}}+\frac{3 e}{4 (b d-a e)^2 (a+b x)^2 (d+e x)^{3/2}}-\frac{21 e^2}{8 (b d-a e)^3 (a+b x) (d+e x)^{3/2}}-\frac{105 b e^3}{8 (b d-a e)^5 \sqrt{d+e x}}+\frac{105 b^{3/2} e^3 \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 (b d-a e)^{11/2}}\\ \end{align*}

Mathematica [C]  time = 0.0173553, size = 52, normalized size = 0.26 \[ -\frac{2 e^3 \, _2F_1\left (-\frac{3}{2},4;-\frac{1}{2};-\frac{b (d+e x)}{a e-b d}\right )}{3 (d+e x)^{3/2} (a e-b d)^4} \]

Antiderivative was successfully verified.

[In]

Integrate[1/((d + e*x)^(5/2)*(a^2 + 2*a*b*x + b^2*x^2)^2),x]

[Out]

(-2*e^3*Hypergeometric2F1[-3/2, 4, -1/2, -((b*(d + e*x))/(-(b*d) + a*e))])/(3*(-(b*d) + a*e)^4*(d + e*x)^(3/2)
)

________________________________________________________________________________________

Maple [A]  time = 0.212, size = 319, normalized size = 1.6 \begin{align*} -{\frac{2\,{e}^{3}}{3\, \left ( ae-bd \right ) ^{4}} \left ( ex+d \right ) ^{-{\frac{3}{2}}}}+8\,{\frac{{e}^{3}b}{ \left ( ae-bd \right ) ^{5}\sqrt{ex+d}}}+{\frac{41\,{e}^{3}{b}^{4}}{8\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}} \left ( ex+d \right ) ^{{\frac{5}{2}}}}+{\frac{35\,{e}^{4}{b}^{3}a}{3\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}} \left ( ex+d \right ) ^{{\frac{3}{2}}}}-{\frac{35\,{b}^{4}d{e}^{3}}{3\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}} \left ( ex+d \right ) ^{{\frac{3}{2}}}}+{\frac{55\,{e}^{5}{b}^{2}{a}^{2}}{8\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}}\sqrt{ex+d}}-{\frac{55\,{e}^{4}{b}^{3}ad}{4\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}}\sqrt{ex+d}}+{\frac{55\,{e}^{3}{b}^{4}{d}^{2}}{8\, \left ( ae-bd \right ) ^{5} \left ( bxe+ae \right ) ^{3}}\sqrt{ex+d}}+{\frac{105\,{e}^{3}{b}^{2}}{8\, \left ( ae-bd \right ) ^{5}}\arctan \left ({b\sqrt{ex+d}{\frac{1}{\sqrt{ \left ( ae-bd \right ) b}}}} \right ){\frac{1}{\sqrt{ \left ( ae-bd \right ) b}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(e*x+d)^(5/2)/(b^2*x^2+2*a*b*x+a^2)^2,x)

[Out]

-2/3*e^3/(a*e-b*d)^4/(e*x+d)^(3/2)+8*e^3/(a*e-b*d)^5*b/(e*x+d)^(1/2)+41/8*e^3*b^4/(a*e-b*d)^5/(b*e*x+a*e)^3*(e
*x+d)^(5/2)+35/3*e^4*b^3/(a*e-b*d)^5/(b*e*x+a*e)^3*(e*x+d)^(3/2)*a-35/3*e^3*b^4/(a*e-b*d)^5/(b*e*x+a*e)^3*(e*x
+d)^(3/2)*d+55/8*e^5*b^2/(a*e-b*d)^5/(b*e*x+a*e)^3*(e*x+d)^(1/2)*a^2-55/4*e^4*b^3/(a*e-b*d)^5/(b*e*x+a*e)^3*(e
*x+d)^(1/2)*a*d+55/8*e^3*b^4/(a*e-b*d)^5/(b*e*x+a*e)^3*(e*x+d)^(1/2)*d^2+105/8*e^3*b^2/(a*e-b*d)^5/((a*e-b*d)*
b)^(1/2)*arctan(b*(e*x+d)^(1/2)/((a*e-b*d)*b)^(1/2))

________________________________________________________________________________________

Maxima [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: ValueError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*x+d)^(5/2)/(b^2*x^2+2*a*b*x+a^2)^2,x, algorithm="maxima")

[Out]

Exception raised: ValueError

________________________________________________________________________________________

Fricas [B]  time = 2.56335, size = 3717, normalized size = 18.58 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*x+d)^(5/2)/(b^2*x^2+2*a*b*x+a^2)^2,x, algorithm="fricas")

[Out]

[-1/48*(315*(b^4*e^5*x^5 + a^3*b*d^2*e^3 + (2*b^4*d*e^4 + 3*a*b^3*e^5)*x^4 + (b^4*d^2*e^3 + 6*a*b^3*d*e^4 + 3*
a^2*b^2*e^5)*x^3 + (3*a*b^3*d^2*e^3 + 6*a^2*b^2*d*e^4 + a^3*b*e^5)*x^2 + (3*a^2*b^2*d^2*e^3 + 2*a^3*b*d*e^4)*x
)*sqrt(b/(b*d - a*e))*log((b*e*x + 2*b*d - a*e - 2*(b*d - a*e)*sqrt(e*x + d)*sqrt(b/(b*d - a*e)))/(b*x + a)) +
 2*(315*b^4*e^4*x^4 + 8*b^4*d^4 - 50*a*b^3*d^3*e + 165*a^2*b^2*d^2*e^2 + 208*a^3*b*d*e^3 - 16*a^4*e^4 + 420*(b
^4*d*e^3 + 2*a*b^3*e^4)*x^3 + 63*(b^4*d^2*e^2 + 18*a*b^3*d*e^3 + 11*a^2*b^2*e^4)*x^2 - 18*(b^4*d^3*e - 10*a*b^
3*d^2*e^2 - 53*a^2*b^2*d*e^3 - 8*a^3*b*e^4)*x)*sqrt(e*x + d))/(a^3*b^5*d^7 - 5*a^4*b^4*d^6*e + 10*a^5*b^3*d^5*
e^2 - 10*a^6*b^2*d^4*e^3 + 5*a^7*b*d^3*e^4 - a^8*d^2*e^5 + (b^8*d^5*e^2 - 5*a*b^7*d^4*e^3 + 10*a^2*b^6*d^3*e^4
 - 10*a^3*b^5*d^2*e^5 + 5*a^4*b^4*d*e^6 - a^5*b^3*e^7)*x^5 + (2*b^8*d^6*e - 7*a*b^7*d^5*e^2 + 5*a^2*b^6*d^4*e^
3 + 10*a^3*b^5*d^3*e^4 - 20*a^4*b^4*d^2*e^5 + 13*a^5*b^3*d*e^6 - 3*a^6*b^2*e^7)*x^4 + (b^8*d^7 + a*b^7*d^6*e -
 17*a^2*b^6*d^5*e^2 + 35*a^3*b^5*d^4*e^3 - 25*a^4*b^4*d^3*e^4 - a^5*b^3*d^2*e^5 + 9*a^6*b^2*d*e^6 - 3*a^7*b*e^
7)*x^3 + (3*a*b^7*d^7 - 9*a^2*b^6*d^6*e + a^3*b^5*d^5*e^2 + 25*a^4*b^4*d^4*e^3 - 35*a^5*b^3*d^3*e^4 + 17*a^6*b
^2*d^2*e^5 - a^7*b*d*e^6 - a^8*e^7)*x^2 + (3*a^2*b^6*d^7 - 13*a^3*b^5*d^6*e + 20*a^4*b^4*d^5*e^2 - 10*a^5*b^3*
d^4*e^3 - 5*a^6*b^2*d^3*e^4 + 7*a^7*b*d^2*e^5 - 2*a^8*d*e^6)*x), 1/24*(315*(b^4*e^5*x^5 + a^3*b*d^2*e^3 + (2*b
^4*d*e^4 + 3*a*b^3*e^5)*x^4 + (b^4*d^2*e^3 + 6*a*b^3*d*e^4 + 3*a^2*b^2*e^5)*x^3 + (3*a*b^3*d^2*e^3 + 6*a^2*b^2
*d*e^4 + a^3*b*e^5)*x^2 + (3*a^2*b^2*d^2*e^3 + 2*a^3*b*d*e^4)*x)*sqrt(-b/(b*d - a*e))*arctan(-(b*d - a*e)*sqrt
(e*x + d)*sqrt(-b/(b*d - a*e))/(b*e*x + b*d)) - (315*b^4*e^4*x^4 + 8*b^4*d^4 - 50*a*b^3*d^3*e + 165*a^2*b^2*d^
2*e^2 + 208*a^3*b*d*e^3 - 16*a^4*e^4 + 420*(b^4*d*e^3 + 2*a*b^3*e^4)*x^3 + 63*(b^4*d^2*e^2 + 18*a*b^3*d*e^3 +
11*a^2*b^2*e^4)*x^2 - 18*(b^4*d^3*e - 10*a*b^3*d^2*e^2 - 53*a^2*b^2*d*e^3 - 8*a^3*b*e^4)*x)*sqrt(e*x + d))/(a^
3*b^5*d^7 - 5*a^4*b^4*d^6*e + 10*a^5*b^3*d^5*e^2 - 10*a^6*b^2*d^4*e^3 + 5*a^7*b*d^3*e^4 - a^8*d^2*e^5 + (b^8*d
^5*e^2 - 5*a*b^7*d^4*e^3 + 10*a^2*b^6*d^3*e^4 - 10*a^3*b^5*d^2*e^5 + 5*a^4*b^4*d*e^6 - a^5*b^3*e^7)*x^5 + (2*b
^8*d^6*e - 7*a*b^7*d^5*e^2 + 5*a^2*b^6*d^4*e^3 + 10*a^3*b^5*d^3*e^4 - 20*a^4*b^4*d^2*e^5 + 13*a^5*b^3*d*e^6 -
3*a^6*b^2*e^7)*x^4 + (b^8*d^7 + a*b^7*d^6*e - 17*a^2*b^6*d^5*e^2 + 35*a^3*b^5*d^4*e^3 - 25*a^4*b^4*d^3*e^4 - a
^5*b^3*d^2*e^5 + 9*a^6*b^2*d*e^6 - 3*a^7*b*e^7)*x^3 + (3*a*b^7*d^7 - 9*a^2*b^6*d^6*e + a^3*b^5*d^5*e^2 + 25*a^
4*b^4*d^4*e^3 - 35*a^5*b^3*d^3*e^4 + 17*a^6*b^2*d^2*e^5 - a^7*b*d*e^6 - a^8*e^7)*x^2 + (3*a^2*b^6*d^7 - 13*a^3
*b^5*d^6*e + 20*a^4*b^4*d^5*e^2 - 10*a^5*b^3*d^4*e^3 - 5*a^6*b^2*d^3*e^4 + 7*a^7*b*d^2*e^5 - 2*a^8*d*e^6)*x)]

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*x+d)**(5/2)/(b**2*x**2+2*a*b*x+a**2)**2,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [B]  time = 1.22311, size = 576, normalized size = 2.88 \begin{align*} -\frac{105 \, b^{2} \arctan \left (\frac{\sqrt{x e + d} b}{\sqrt{-b^{2} d + a b e}}\right ) e^{3}}{8 \,{\left (b^{5} d^{5} - 5 \, a b^{4} d^{4} e + 10 \, a^{2} b^{3} d^{3} e^{2} - 10 \, a^{3} b^{2} d^{2} e^{3} + 5 \, a^{4} b d e^{4} - a^{5} e^{5}\right )} \sqrt{-b^{2} d + a b e}} - \frac{315 \,{\left (x e + d\right )}^{4} b^{4} e^{3} - 840 \,{\left (x e + d\right )}^{3} b^{4} d e^{3} + 693 \,{\left (x e + d\right )}^{2} b^{4} d^{2} e^{3} - 144 \,{\left (x e + d\right )} b^{4} d^{3} e^{3} - 16 \, b^{4} d^{4} e^{3} + 840 \,{\left (x e + d\right )}^{3} a b^{3} e^{4} - 1386 \,{\left (x e + d\right )}^{2} a b^{3} d e^{4} + 432 \,{\left (x e + d\right )} a b^{3} d^{2} e^{4} + 64 \, a b^{3} d^{3} e^{4} + 693 \,{\left (x e + d\right )}^{2} a^{2} b^{2} e^{5} - 432 \,{\left (x e + d\right )} a^{2} b^{2} d e^{5} - 96 \, a^{2} b^{2} d^{2} e^{5} + 144 \,{\left (x e + d\right )} a^{3} b e^{6} + 64 \, a^{3} b d e^{6} - 16 \, a^{4} e^{7}}{24 \,{\left (b^{5} d^{5} - 5 \, a b^{4} d^{4} e + 10 \, a^{2} b^{3} d^{3} e^{2} - 10 \, a^{3} b^{2} d^{2} e^{3} + 5 \, a^{4} b d e^{4} - a^{5} e^{5}\right )}{\left ({\left (x e + d\right )}^{\frac{3}{2}} b - \sqrt{x e + d} b d + \sqrt{x e + d} a e\right )}^{3}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*x+d)^(5/2)/(b^2*x^2+2*a*b*x+a^2)^2,x, algorithm="giac")

[Out]

-105/8*b^2*arctan(sqrt(x*e + d)*b/sqrt(-b^2*d + a*b*e))*e^3/((b^5*d^5 - 5*a*b^4*d^4*e + 10*a^2*b^3*d^3*e^2 - 1
0*a^3*b^2*d^2*e^3 + 5*a^4*b*d*e^4 - a^5*e^5)*sqrt(-b^2*d + a*b*e)) - 1/24*(315*(x*e + d)^4*b^4*e^3 - 840*(x*e
+ d)^3*b^4*d*e^3 + 693*(x*e + d)^2*b^4*d^2*e^3 - 144*(x*e + d)*b^4*d^3*e^3 - 16*b^4*d^4*e^3 + 840*(x*e + d)^3*
a*b^3*e^4 - 1386*(x*e + d)^2*a*b^3*d*e^4 + 432*(x*e + d)*a*b^3*d^2*e^4 + 64*a*b^3*d^3*e^4 + 693*(x*e + d)^2*a^
2*b^2*e^5 - 432*(x*e + d)*a^2*b^2*d*e^5 - 96*a^2*b^2*d^2*e^5 + 144*(x*e + d)*a^3*b*e^6 + 64*a^3*b*d*e^6 - 16*a
^4*e^7)/((b^5*d^5 - 5*a*b^4*d^4*e + 10*a^2*b^3*d^3*e^2 - 10*a^3*b^2*d^2*e^3 + 5*a^4*b*d*e^4 - a^5*e^5)*((x*e +
 d)^(3/2)*b - sqrt(x*e + d)*b*d + sqrt(x*e + d)*a*e)^3)